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ABSTRACT 

The coverage and quality of remotely sensed upper-tropospheric moisture parameters have improved considerably 
with the deployment of a new generation of operational geostationary meteorological satellites: GOES-8/9 and GMS-5. 
The GOES-8/9 water vapor imaging capabilities have increased as a result of improved radiometric sensitivity and higher 
spatial resolution. The addition of a water vapor sensing channel on the latest GMS permits nearly global viewing of 
upper-tropospheric water vapor (when joined with GOES and Meteosat) and enhances the commonality of geostation- 
ary meteorological satellite observing capabilities. Upper-tropospheric motions derived from sequential water vapor 
imagery provided by these satellites can be objectively extracted by automated techniques. Wind fields can be deduced 
in both cloudy and cloud-free environments. In addition to the spatially coherent nature of these vector fields, the 
GOES-8/9 multispectral water vapor sensing capabilities allow for determination of wind fields over multiple tropo- 
spheric layers in cloud-free environments. This article provides an update on the latest efforts to extract water vapor 
motion displacements over meteorological scales ranging from subsynoptic to global. The potential applications of these 
data to impact operations, numerical assimilation and prediction, and research studies are discussed. 


1. Introduction 

Two of the most important challenges facing the at- 
mospheric remote sensing community today are 1 ) trans- 
forming information-wealthy imagery into quantitative 
parameters and 2) ensuring these parameters advance our 
understanding and prediction of meteorological phenom- 
ena beyond existing methods of observing the atmo- 
sphere. The latter point is of particular significance 
given the advances in numerical weather prediction 
(NWP) over the past decade. Further improvements 
in prediction will partially rely on superior measure- 
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ments of the initial atmospheric state and optimized 
assimilation of these measurements. This is of particu- 
lar relevance to remotely sensed observations when 
considering regions of the globe that are traditionally 
void of conventional meteorological data. 

Water vapor radiance measurements from geosta- 
tionary satellites have been used extensively in a quali- 
tative mode by analysts and forecasters (Weldon and 
Holmes 1991). Early attempts to quantify the data by 
employing manual methods to track moisture features 
in animated image sequences met with modest suc- 
cess (Hayden and Stewart 1987). More recent stud- 
ies have demonstrated the ability of automated 
extraction techniques to produce spatially coherent 
water vapor wind vector (WVWV) fields comparable 
in quality to operational cloud-tracked winds (Laurent 
1993; Velden 1993; Holmlund 1993). 

Development of automated wind-extraction algo- 
rithms at the University of Wisconsin — Madison Co- 
operative Institute for Meteorological Satellite Studies 
(CIMSS) has progressed to the stage of operational 
implementation (Nieman et al. 1997). In collaboration 
with the National Oceanic and Atmospheric Admin- 
istration (NOAA) National Environmental Satellite 
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Data and Information Service (NESDIS), efforts at 
CIMSS have focused on the research, development, 
and demonstration of these algorithms toward the goal 
of routine inclusion of the WVWV measurements into 
the operational data stream and numerical weather 
prediction. Attributable to the advances in computing 
power, algorithm automation, and communications, 
high-density WVWV fields from multiple geostation- 
ary satellites each covering near full-disk viewing are 
achievable in a time frame and frequency commen- 
surate with operational requirements. 

This article describes the general characteristics of 
satellite-derived WVWV, the methodology employed 
at CIMSS to derive the vector fields, and the currently 
attainable accuracies. Several examples of applica- 
tions will be presented, as well as initial numerical 
model assimilation impact results. This is followed by 
a brief discussion of the potential benefits to the ob- 
serving system, the limitations and challenges in ef- 
fective use of the data, and future prospects. 

2. Background 

The quest for quantitative wind information over 
oceanic regions lacking in conventional observations 
was the prime motivation for the development of 
methods to track clouds in sequential satellite imag- 
ery. These “cloud-drift winds” became an important 
component of the global observing system during the 
early 1980s. However, an obvious limitation of this 
product resulted from the fact that vast cloud-free re- 
gions were left unsampled. About this time, 6-7-//m 
water vapor absorption band channels became avail- 
able on geostationary satellites (Meteosat and GOES). 
Imagery from these observations revealed atmo- 
spheric motion even in cloud-free areas, and research 
with these data led Eigen willig and Fischer (1982) to 
conclude that feature tracking was possible within a 
reasonable accuracy. This was supported in a study 
by Stewart et al. (1985) using GOES imagery and a 
state-of-the-art tracking system (McIDAS, Man- 
computer Interactive Data Access System). 

These early attempts at WVWV employed basic tech- 
niques similar to cloud tracking, in which a human op- 
erator selected a target feature in an image and sub- 
jectively followed it through subsequent imagery to 
determine a displacement. The heights were crudely as- 
signed utilizing basic radiative transfer formulae to trans- 
form water vapor (WV) pixel radiance to a brightness 
temperature, which was then matched to a collocated 


temperature profile (model guess or climatology). 
While WVWV fields were successfully demonstrated 
(LeMarshall et al. 1985; Hayden and Stewart 1987), the 
subjectivity and labor-intensive nature of the technique 
prevented the product from becoming fully operational. 

In the late 1980s and early 1990s, advancement in 
the automation of satellite feature tracking tech- 
niques (Merrill et al. 1991) was partially successful 
toward easing the production constraints on WVWV 
fields. However, the quality control (QC) of the 
WVWV fields still demanded intensive manual inter- 
vention, which restricted operational utility. In addi- 
tion, numerical model forecast accuracies were 
increasing rapidly, and the quality of the GOES-7 
WVWV observations was insufficient to make signifi- 
cant impacts (Velden 1990). By contrast, the European 
community (European Space Operations Center) de- 
veloped an automated WVWV extraction algorithm 
using higher quality Meteosat observations. This prod- 
uct became operational in 1993, although the QC is 
very strict (Laurent 1993; Holmlund 1995). 

Recognizing the need for improved imaging and 
sounding capabilities from geostationary orbit, NOAA 
began development of its next generation of satellites, 
GOES I-M, in 1985 (Menzel and Purdom 1994). The 
first two of these GOES were successfully deployed 
in 1994 and 1995. This new generation features ad- 
vanced sensors enabling improved spatial resolution and 
signal-to-noise ratios. Specifically, the water vapor 
channels (imager and sounder) are a significant up- 
grade from the previous GOES (Table 1 ). A qualitative 
comparison of GOES-7 and GOES -8 water vapor im- 
agery is shown in Fig. 1 and illustrates the signal-to-noise 
improvements associated with the new generation. 

The remainder of this article will focus on the cur- 
rent ability to extract multispectral (imager and sounder) 
WVWV information from geostationary satellites. It will 
be shown that the combination of improved radiomet- 
ric observations, advanced and automated wind extrac- 
tion algorithms, and an optimized processing strategy 
can lead to a product that has the potential to make 
an important contribution to the upper-tropospheric 
observing system. The WVWV applications cover 
meteorological scales ranging from meso- to climate 
and can benefit both research and operations. 

3. Methodology 

The procedure for deriving WVWV at CIMSS par- 
allels the basic technique employed for cloud track- 
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Fig. 1. water vapor imagery over the south-central United States trom 

GOES-7 (left) and GOES-8 (right) at 1 300 UTC 3 1 May 1 994. The color enhancement 
reflects the brightness temperature range with stronger yellows indicating drier air. 


ing. The fully automated algorithm is 
housed within McIDAS-X (UNIX) 
architecture. A description of the sys- 
tem can be found in Merrill (1989) 
and Merrill et al. (1991), with the lat- 
est upgrades (including specialized 
image registration) in Nieman et al. 
(1997). Only a brief account will be 
given here, with attention given to 
features unique to the derivation of 

wvwv. 


a. Tracer selection and tracking 

metric 

The traditional selection process 
for targeting features in the case of 
tracking clouds is based on the isola- 
tion of the lowest (cloud) brightness temperatures 
(BT) within a pixel array (target box) of a size selected 
by the operator. Local bidirectional gradients are com- 
puted and compared with empirically determined 
thresholds to identify those features exhibiting suffi- 
cient variability to be classified as targets (Merrill 
1989). 

Water vapor structure in cloud-free environments 
demands a different strategy. Features exhibiting the 
strongest gradient profiles may not be confined to the 
coldest BT, as in clouds. In the CIMSS algorithm, 
targets are identified by evaluating the bidirectional 
gradients surrounding each pixel in the target array 
and selecting the maximum value (exceeding thresh- 
olds). The operator can control the density of targets 
by selecting the size of the target box and setting the 
gradient thresholds. This flexibility is important when 
considering different meteorological scales to be mea- 
sured and for computational efficiency. 

In practice, WVWV targets are selected in both 
cloudy and cloud-free regions. Experience has shown 
that the density of higher-quality targets is greater in 
cloudy and/or very moist regimes exhibiting high con- 
trast. This is reflected in the resulting WVWV wind 
field characteristics, which are discussed in the next 
section. 

The tracking metric currently employs a simple 
search for the minimum in the sum of squares of the 
radiance differences between the target and search 
arrays in subsequent images (usually two subsequent 
images at 30-min intervals). Fully cross-correlational 
algorithms are more computationally demanding but 
are being investigated. A model guess forecast of the 
upper-level wind is employed to economize on the 


search areas. In practice, two subsequent images are 
utilized with the target image to derive two displace- 
ment vectors. If the vectors survive consistency 
checks, an average is calculated and becomes the rep- 


Table 1. Performance specifications and spatial resolution of 
current geostationary satellite water vapor channels. NEDT is the 
noise equivalent temperature difference (K) evaluated at 230 K. 
Numbers in parentheses indicate central wavelengths. 


NEDT 

Resolution (km) 

Imagers 

GOES-9 (6.7) 

0.11 

8 

GOES-8 (6.7) 

0.22 

8 

GOES -7(6.7) 

1.00 

14 

Meteosat-5 (6.5) 

0.40 

5 

GMS-5 (6.9) 

0.40 

5 

Sounders 

GOES-9 (7.3) 

0.33 

8 

GOES-9 (7.0) 

0.27 

8 

GOES-8 (7.3) 

0.51 

8 

GOES-8 (7.0) 

0.40 

8 

GOES-7 (7.0) 

3.70 

14 
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resentative WVWV. Experimentation has shown that 
30-min interval imagery is slightly superior to hourly 
intervals for optimal WVWV tracking. 

b. Height assignment and quality control 

A long-standing concern regarding WVWV is the 
height assignment. Particularly in cloud-free regions, 
the radiometric signal from pure water vapor struc- 
tures is a result of emittance over a finite layer. This 
can be complicated further by radiance contributions 
from multiple moist layers (Weldon and Holmes 
1991). The challenge is to assign a height level that 
best represents the motion of the feature being tracked. 

The strategy employed in the CIMSS algorithm for 
height assignment is a two-step approach. The first is 
the conversion of instrument measured radiances into 
BT. This requires an accurate definition of the spec- 
tral response of the instrument. Figure 2 shows the 
spectral response functions for GOES-8, Meteosat-5 , 
and GMS-5 water vapor channels. It is important to 
note that the different instruments are sensing radia- 
tion from slightly different parts of the H.O absorp- 
tion band. Sampling the center of the absorption band 
(around 6.3 ja m) yields radiation from the upper lev- 
els of the troposphere (i.e., radiation from below has 
already been absorbed by the water vapor). Sampling 
away from the center of the absorption band yields 
radiation from lower levels of the troposphere. 
Figure 3 shows the weighting functions (which denote 
energy contribution profiles) for the GOES-8 imager 
and sounder H,0 sensitive bands, as well as the 
Meteosat-5 and GMS-5 water vapor bands. The 



Fig. 2. Water vapor channel spectral response functions for 
GOES-8, Meteosat-5 , and GMS-5. The scale at the left indicates 
a normalized measure of the relative energy contribution per 
spectral region. 


GOES-8 sounder 7.3-jum band senses deeper into the 
atmosphere than the other bands, as it is farther away 
from the H.O absorption band center. The GMS-5 
6.9-^im band senses a similar layer of the atmosphere 
as the GOES-8 sounder 7.0-/tm band. 

Once a target has been selected for tracking, the BT 
in the target area are averaged and this value is 
matched with a collocated model guess temperature 
profile. The level of optimum fit is used to assign an 
initial pressure height. As indicated in Fig. 3, these 
heights are typically within the range of 250-600 mb. 

The second step of the height assignment procedure 
involves an important and unique component of the 
CIMSS winds processing system and is automatically 




Fig. 3. Weighting functions associated with geostationary 
satellite water vapor channels in a typical subtropical atmosphere: 
GOES-8 (top), Meteosat-5 and GMS-5 (bottom). Numbers in 
parentheses refer to central wavelengths. Peaks in the functions 
correspond to tropospheric layers of greatest energy contribution. 
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invoked following the completed computation of the 
WVWV field. The initial WVWV heights are reevalu- 
ated as part of the postprocessing algorithm (objec- 
tive QC procedure outlined below). The basic 
rationale for this height reassessment is to seek the as- 
signment of a level that best fits other multivariate 
information from conventional data, neighboring 
WVWV and cloud-tracked winds, numerical model 
output, etc. In this sense, the procedure has an element 
of data assimilation (Hayden and Velden 1991). 

The postprocessing algorithm involves a two-stage, 
three-dimensional objective analysis (Hayden and 
Purser 1995) of the wind field using background in- 
formation from a numerical model forecast. The 
procedure is schematically shown in Fig. 4. The ob- 
jective analysis employed is an adaptation of the re- 
cursive filter (RF) analysis described by Hayden and 
Purser (1988). It is a basic successive approximation 
method but contains a unique feature of locally vary- 
ing scaling, which gives it greater flexibility over 
nonhomogeneous data. The general properties of the 
RF, and its relationship to a characteristic scale, are 
given in Hayden and Purser (1995). 

The first stage of the postprocessing algorithm pro- 
vides a preliminary wind analysis using the WVWV 
at their initially assigned pressures, any ancillary re- 
ports (i.e., cloud-drift winds), and pseudodata from a 
numerical model 6- or 12-h forecast. At this point, an 
adjustment to the WVWV speed is applied based on 
the well-known slow bias in upper-tropospheric 
cloud-drift winds. It is suspected that this bias is a 
result of clouds and water vapor motions not being a 
perfect indicator of the instantaneous wind, especially 
in higher wind regimes. This slow bias is mitigated 
by an empirically determined correction term (partly 
based on statistics from the European Centre for 
Medium-Range Weather Forecasts), which is applied 
to WVWV with speeds greater than 10 m s _1 . The ad- 
justment increments each vector by 8% of vector 
speed. 

The WVWV heights are then adjusted by conduct- 
ing a vertical search aimed at minimizing a simple 
variational penalty function evaluated using the ini- 
tial wind analyses and the model forecast: 
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where V = velocity, T = temperature, P = pressure, 
dd = direction, and s = speed. Here, subscript m re- 
fers to a WVWV measurement, / and j are horizontal 
dimensions in the analysis, and k is the vertical level. 
The denominators F x are weights defined by the op- 
erator that may be varied to emphasize the vector ve- 
locity, temperature (target BT), pressure, direction, 
and speed terms in the penalty function. The vertical 
search for a best fit (and hence, the maximum allow- 
able reassignment) is limited to ±150 mb from the 
initial assignment. The height reassignment may fail 
if no minimum is found or if the minimum exceeds 
empirically determined thresholds. In these instances 
the vector is rejected. 

The full RF analysis is repeated in stage two of the 
postprocessing algorithm using the reassigned and 
accelerated WVWV and output of the first analysis 
as a background field. The output of this analysis pro- 
vides quality estimates (flags) for each vector based 
on the local quality of the analysis and the fit of the 
observation to that analysis. Vectors that do not at- 
tain an empirically determined threshold quality are 
rejected. 

In regard to the threshold, a trade-off exists in the 
choice of the maximum allowable difference between 
the background field (which may be highly influenced 
by the model guess winds) and the vectors: if only a 
small difference is allowed, potential vector informa- 
tion is reduced; if a large difference is allowed, po- 
tential information content is increased but so are 
potential errors. The threshold currently employed in 
the operational NESDIS algorithm (a version of the 
CIMSS algorithm) is rather conservative and is based 
on an error analysis against collocated rawinsonde 
information. However, this threshold is not optimal 
in meteorological situations subject to difficult model 
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Fig. 4. Schematic of the objective quality control procedure 
employed in the CIMSS automated water vapor wind vector 
(WVWV) processing algorithm. 
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forecasts, such as hurricanes or strong extratropical 
events. In this regard, a wide variety of options are 
available for regulating the analysis (such as the back- 
ground field weights), the penalty function, and the 
final quality estimates. For WVWV, these options 
have been empirically tuned over several years of 
experimenting and application. 

The precision attainable from current geostation- 
ary satellite water vapor channels allows WVWV to 
be processed at high density, providing a uniform, 
spatially coherent field of coverage. With GOES, 
WVWV from two sounder channels can be added 
to the analysis, resulting in augmentation of three- 
dimensional coverage (see next section). To take full 
advantage of this coverage, the postprocessing algo- 
rithm is tunable in the following ways: 1) the objec- 
tive analysis horizontal resolution (increment) can be 
increased to capture finer-scale features, 2) the model 
background field can be downweighted in the analy- 
sis, 3) a tighter fit of each WVWV to the analysis can 
be demanded due to the abundance of “neighbor 
checking,” 4) added weight can be given to the ve- 
locity term in the penalty function, 5) vertical coupling 
requirements in the three-dimensional analysis can be 
made less strict, and 6) the rejection criteria based on 
the final quality flags can be relaxed. The final point 
indicates confidence in the ability of the WVWV 
fields to capture and retain details not present in the 
model forecast background field. 

The result of this postprocessing strategy is an 
objectively based, quality-controlled WVWV field 
with vector heights assigned to pressure levels based 
on both radiometric properties and optimum fit with 
ancillary information. While in nature it is rare that a 
particular single level will perfectly represent a layer, 
our analysis of WVWV height assignment quality 
compared with collocated rawinsondes generally in- 
dicates good agreement and placement confidence 
within 50 mb in most cases. WVWV accuracies and 
verification statistics are presented in section 5. 

4. WVWV characteristics 

An example of the attainable coverage from 
GOES-8 6.7-jum WVWV is shown in Fig. 5. This field 
was produced using the CIMSS algorithm on a desk- 
top workstation in about 1 0 min. The coverage domain 
in this example represents a small portion of what is 
possible from geostationary viewing. Featured is 
Tropical Storm Pablo and its environment. The 


WVWV have been processed at high density, using 
the tuned QC strategy described above. Several prop- 
erties of the WVWV product are apparent from this 
example: 1) the coverage is uniform and spatially co- 
herent, 2) the vertical distribution of information is 
spread over the 150-500-mb portion of the tropo- 
sphere, 3) large-scale and synoptic features are well 
represented, and 4) smaller-scale features such as the 
storm outflow and adjacent environmental features are 
captured. 

The last point is of particular importance, as this 
scale is often undersampled by the existing observa- 
tional network and as a result is frequently not well 
represented in global objective analyses. For example, 
the narrow band of high-level 50-60-kt westerly 
winds just upstream (west) of Pablo was not captured 
by any operationally available data and not analyzed 
by the National Centers for Environmental Prediction 
(NCEP) operational global model. Pablo was mov- 
ing steadily to the west and forecasted to slowly in- 
tensify. However, as Pablo encountered these winds 
over the next 24 h, the strong shear disrupted the 
storm circulation and Pablo was subsequently down- 
graded to a depression. This illustrates how the 
WVWV can be utilized in a qualitative mode, since 
knowledge of the impending hostile upper-level en- 
vironment could have played a role in the intensity 
forecast. 

The vertical distribution of WVWV is dependent 
on several key factors: 1) characteristics of the spec- 
tral channel, such as central wavelength and band 
width; 2) airmass characteristics (including the 
amount of cloudiness and vertical profile of moisture); 
3) season (less moisture in the cool season leads to 
lower mean heights); and 4) latitude (higher moisture 
content and tropopause in the Tropics leads to higher 
mean heights). In regard to the first point, the current 
operational geostationary weather satellites with wa- 
ter vapor channels (GOES, Meteosat, and GMS) have 
similar traits (Fig. 3 and Table 1 ), which is advanta- 
geous for global data assimilation and future climate 
research studies. 

In addition to the 6.7-jum channel on the imager, 
GOES-819 house an atmospheric sounder that in- 
cludes two water vapor absorption channels centered 
at 7.0 and 7.3 fm\ (Menzel and Purdom 1994). These 
sounder channels can be employed as surrogate im- 
agers to track water vapor features radiating from 
lower layers of the troposphere (Fig. 3). An example 
of the WVWV coverage attainable from these two 
channels ( GOES-8 ) is illustrated in Fig. 6. Of note is 
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the spatial coherency in the WVWV fields derived however, a programmable scanning strategy allows 
from each channel. A histogram analysis of WVWV the sounder to view targeted areas of meteorological 
heights for this case is shown in Fig. 7. The heights interest. 

are partitioned into 1 2.5-mb bins and graphed to show The coverage attainable from WVWV is quite uni- 
the vertical distribution of WVWV for the 6.7- (im- form in most cases, since the observed displacement 
ager channel 3), 7.0- (sounder channel 1 1 ), and 7.3- of a large portion of the targeted features is a result 
(sounder channel 10) /urn channels. Quasi-independent of pure moisture advection and adequately represents 
peaks in information are found in discrete layers; near the instantaneous wind (after the bias correction). The 
250 mb (cloud), 300-350 mb (channel 3), 350- coverage is especially complete in the summertime 
450 mb (channel 1 1 ), and 450-550 mb (channel 1 0). Tropics due to the presence of abundant moisture and 
While there is some redundancy in terms of cover- other characteristics mentioned previously. There are, 
age in the vertical, additional information is gained however, situations when extraction of useful WVWV 
from the sounder channels, especially in the mid- information is difficult, if not impossible. These situ- 
tropospheric cloud-free regions. Integration of the ations are summarized as follows: 1 ) very dry air 
sounder WVWV is important to the objective, masses lacking in moisture features, which prohibits 
three-dimensional analysis, which is part of the in- good target selection; 2) regions ot strong vertical 
ternal quality control discussed in the previous sec- motion, in which a coherent radiance signal can be 
tion. The sounder coverage domain is restricted due significantly attenuated over the tracking interval; 
to observing constraints (Menzel and Purdom 1 994); 3) uniform moisture fields lacking in contrast (flat gra- 



Fig. 5. High-density GOES-8 water vapor winds (6.1 -/urn channel) during Tropical Storm Pablo (lower right) on 7 October 1995. 
Vectors are in kt and colors refer to assigned heights (mb) with the legend in the upper right-hand comer. 
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and instrument characteristics 
such as detector precision and 
resolution. The next section will 
evaluate the accuracy of WVWV 
generated by the CIMSS track- 
ing system through comparisons 
with rawinsonde observations. 


5. WVWV verification 

The evaluation of WVWV 
must take into account both 
qualitative and quantitative 
measures. Certainly, to observa- 
tionalists and operational fore- 
casters, the complementary 
nature of the WVWV to cloud- 
drift winds in cloudy regions is 
an important enhancement to the 
qualitative analysis of the upper- 
tropospheric wind field. For ex- 
ample, it was shown in the last 
section how the WVWV pro- 
vided information on the envi- 
ronment of Tropical Storm 
Pablo that may have been use- 
ful in the intensity forecast. 
Another good example is the 
WVWV fields being derived 
routinely from GMS-5 in a dem- 
Fig. 6. Example of wind coverage that can be achieved from tracking features in the onstration mode by CIMSS for 
7.0- and 7.3-/im channels from the GOES-8 sounder. Yellow (beige) vectors (kt) indicate the Joint Typhoon Warning 
winds with assigned heights above (below) 400 mb. The winds are plotted over their Center (JTWC) at Guam These 

K data are being integrated into the 

JTWC subjective upper-level 
analyses, and feedback has been 

dients); 4) strong vertical shear zones, which can de- very positive (further details in section 7). 
form targets over the tracking interval; and 5) rapidly Quantitative assessment of the WVWV product is 
amplifying extratropical troughs, which contain a possible from statistical analyses and impact on NWP. 
large phase speed component (over the tracking in- The latter will be addressed in the next section. A tra- 
terval, the observed WVWV displacement is not a ditional method of validation is matching observations 
good representation of the instantaneous flow), with collocated rawinsondes. These statistics can pro- 
Fortunately, the above conditions represent relatively vide a fixed measure of product quality over time and 
rare situations when considering full-disk (roughly 65° can be employed in determining observation weight 
radius from satellite subpoint) coverage, and experi- in objective data assimilation, 
ence shows good coverage is attainable over the large Validation statistics for WVWV from the GOES-8 
majority of meteorological conditions. 6.7-jtm channel are shown in Table 2. These WVWV 

In summary, the coverage and quality of WVWV were produced by NESDIS using the CIMSS algo- 
will vary depending on atmospheric factors such rithm during the months of October-December 1995 
as moisture availability and structure, vertical shear, [at the time this article was submitted, the NESDIS 
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Fig. 7. Typical distribution of water vapor wind vector heights from a 
multispectral GOES-8 dataset: blue = channel 3 (6.7 ji m), red = channel 1 1 
(7.0 jm\), and green = channel 10 (7.3 Jim). The frequency is the number of winds 
(x 10 2 ) per 12.5-mb altitude bin. 


WVWV were being produced in a quasi- 
operational mode and evaluated by 
NCEP-EMC (Environmental Modeling 
Center) for consideration as input to 
NCEP operational models]. Table 2 
shows that the overall quality controlled 
WVWV root-mean-square errors (rmse) 
during this period were 7.47 m s _l , with 
a bias of -0.26 m s' 1 . These rmse inher- 
ently include errors in rawinsonde mea- 
surement and reporting, which are 
typically around 3 m s _l (Hoehne 1980). 

The statistics are further stratified by 
vector height and location. Generally, the 
best results are achieved in the upper tro- 
posphere (above 400 mb) and in the 
Tropics, for reasons related to character- 
istics discussed in the previous section. 

It is insightful to compare the statis- 
tics of WVWV with operational cloud 
drift winds from GOES-8 produced by 
NESDIS over the same period (Table 3). 

Only the high- and midlevel stratifica- 
tions can be evaluated, since there are no 
low-level (below 700 mb) WVWV. 
Comparatively, the quality of high-level WVWV 
match well with the operational cloud drift winds 
but is slightly inferior at midlevels (below 400 mb). 
However, the coverage and information in the 
midlevels is almost twice that of the cloud-drift winds. 

The statistics presented in Table 2 are a reflection 
of WVWV produced with a rather stringent internal 
quality control that is being applied in the postpro- 
cessing stage. NESDIS and operational considerations 
demand this approach. However, research at CIMSS is 
suggesting an alternate processing strategy aimed at 
relaxing the QC constraints slightly in favor of retain- 
ing information in key meteorological areas (discussed 
in section 3). The result of this approach yields vec- 
tor rmse that are slightly higher. However, the trade- 
off can be dramatic in terms of retaining vector 
coverage and information in regions of model guess 
deficiencies. For example, as shown earlier, the ap- 
plication of WVWV data to the analyses of tropical 
cyclone wind fields can yield significant information 
on the upper-level environmental circulations affect- 
ing storm motion ( Velden 1 996; Franklin et al. 1 996). 

For purposes of data assimilation, the WVWV ac- 
curacies are related to vector wind speed in Fig. 8. 
Generally speaking, vector rmse increase steadily with 
wind speed. Part of this may be attributed to increased 


errors in rawinsonde measurements at higher wind 
speeds (SRCC 1981). However, if the WVWV rmse 
values are considered as a percent of vector speed, the 
accuracy actually increases in higher wind regimes. 

It was mentioned in the previous section that 
WVWV are also possible from GOES sounder chan- 
nels measuring radiances at central wavelengths of 7.0 
and 7.3 jum, which allow sensing in cloud-free regions 
in the midtroposphere (Fig. 3). Statistical evaluations 
done at CIMSS show these vector rmse to be compa- 
rable to the 6.7-jUm WVWV. Figures 9 and 10 illus- 
trate the vector rmse relative to wind speed for the 7.0- 
and 7.3-jum WVWV, respectively. The magnitude and 
characteristic increase in rmse with respect to vector 
speed is similar to that found with the 6.7-jUm WVWV. 
Comparisons in higher wind speed regimes are more 
limited due to the nature of the sounder WVWV sens- 
ing at relatively lower tropospheric levels (and wind 
speeds). In the midlevels, the lack of good moisture 
signatures or multiple moist layers within the energy 
contribution function layer is more common and poses 
a greater challenge to the retrieval of accurate tracers 
and attendant height assignments. Preliminary statis- 
tical analyses and NWP results from NCEP-EMC are 
indicating a 1-2 m s" 1 slow speed bias in the 7.3-jum 
WVWV in the lowest levels. This is being investigated. 


Bu/fetin of the American Meteorological Society 


181 



12.0 


Table 2. Rawinsonde verification statistics for GOES-8 
6.7-jUm water vapor winds during October-December 1995. The 
collocation match radius does not exceed 2 degrees. Root-mean- 
square vector error (rmse), standard deviation of the vector 
difference (stdv), mean speed difference (bias), rawinsonde speed 
(spd) (all m s -l )> and number of collocated matches (num) are 
shown. High represents WVWV between 100 and 399 mb, and 
mid is 400-699 mb. 



All 

regions 

NH 

extratopics 

Tropics 

SH 

extratropics 

All levels 

Rmse 

7.47 

7.97 

5.61 

8.02 

Stdv 

4.45 

4.71 

3.08 

4.77 

Bias 

-0.26 

-0.32 

-0.16 

0.78 

Spd 

23.04 

26.87 

10.69 

18.97 

Num 

9337 

7030 

2118 

194 

High 

Rmse 

7.45 

7.95 

5.62 

8.37 

Stdv 

4.43 

4.71 

3.09 

4.92 

Bias 

-0.13 

-0.14 

-0.17 

0.74 

spd 

23 38 

27.70 

10.75 

i9 72 

Num 

8378 

6151 

2067 

165 

Mid 

Rmse 

7.93 

8.12 

5.57 

5.69 

Stdv 

4.60 

4.70 

2.63 

3.32 

Bias 

-1.41 

-1.58 

0.18 

0.97 

Spd 

20.15 

21.02 

8.15 

14.66 

Num 

959 

879 

51 

29 


CIMSS is routinely producing WVWV twice a day 
from GMS-5 for analysis by the JTWC. The process- 
ing strategy is geared to a high-density product cov- 
ering the western North Pacific Ocean, with adjusted 
QC constraints as outlined in section 3. An initial 
evaluation of these observations versus rawinsonde 
data is yielding vector quality on par with what is 
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Fig. 8 . (Top): Accuracy (vector rms difference) of b.l-jim 
GOES-8 water vapor wind vectors produced quasi-operationally by 
NESDIS compared to collocated rawinsondes over the United States 
with respect to vector wind speed (3 m s _l bins). (Bottom): Number 
of wind-raw in sonde matches per wind speed bin. The matches 
were collected during the period of October-December 1995. 


being observed with GOES WVWV being produced 
at CIMSS: rmse of 7.0 m s _l , with a bias of-0.5 m s" 1 . 
The GMS-5 WVWV accuracies relative to vector 
speed are shown in Fig. 1 1 . The functional relation- 
ship is quite similar to that for GOES-8 illustrated in 
Figs. 8-10. The Japanese Meteorological Agency 
(JMA) is also experimenting with automated tech- 
niques to produce WVWV from the GMS-5. At the 
time of writing, this product is being evaluated for 
operational use at the JMA. 


6.NWP data sensitivity experiments 

The importance of observations is often measured 
in terms of the impact on objective analyses and nu- 
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Fig. 9. Same as Fig. 8 except for WVWV produced by CIMSS 
from the 7.0-jt/m channel. 


merical prognoses. Previous studies have provided a 
glimpse into the impact of WVWV data on NWP. 
Despite a nonoptimized assimilation process, the data 
were found to have a positive impact on regional fore- 
casts utilizing a limited-area prediction system 
(Velden 1996). However, the impact was modest. The 
WVWV were also assimilated into a low-resolution 
(T62) version of the NCEP global spectral model for 
a period of a month in 1 994. The results were not sta- 
tistically significant, but virtually all parameters ex- 
amined exhibited a small positive impact (Velden 
1996). The strongest positive impact was noted in the 
Tropics, in particular, the rotational component of the 
wind. 

It was shown in Velden et al. (1992) that the in- 
clusion of high-density satellite-derived wind in- 
formation into the analyses of tropical cyclone 
environmental wind fields can effectively reduce the 
error of objective track forecasts. This particular study 


Tablh 3. Rawinsonde verification statistics for GOES-8 IR 
cloud-drift winds during October-December 1995. The colloca- 
tion match radius does not exceed 2 degrees. Root-mean-square 
vector error (rmse), standard deviation of the vector difference 
(stdv), mean speed difference (bias), rawinsonde speed (spd) (all 
m s~‘X and number of collocated matches (num) are shown. High 
represents winds between 100 and 399 mb, mid is 400-699 mb, 
and low is 700-950 mb. 


All 

NH 


SH 


regions 

extratopics 

Tropics i 

extratropics 

All levels 

Rmse 

7.07 

7.50 

5.47 

7.54 

Stdv 

4.26 

4.54 

2.98 

4.85 

Bias 

-0.40 

-0.41 

-0.46 

0.36 

Spd 

19.20 

22.31 

10.48 

16.46 

Num 

2823 

2039 

703 

81 

High 

Rmse 

7.44 

7.88 

5.93 

8.32 

Stdv 

4.41 

4.72 

3.19 

5.12 

Bias 

-0.43 

-0.26 

-0.97 

0.54 

Spd 

22.59 

26.81 

11.81 

20.80 

Num 

1917 

1346 

517 

54 

Mid 

Rmse 

6.34 

6.64 

4.40 

4.22 

Stdv 

3.72 

3.90 

2.04 

2.40 

Bias 

-0.52 

-0.70 

0.61 

-0.70 

Spd 

13.94 

15.18 

6.43 

12.18 

Num 

554 

469 

75 

10 

Low 

Rmse 

5.97 

6.84 

3.33 

6.50 

Stdv 

3.95 

4.52 

L46 

4.56 

Bias 

-0.04 

-0.66 

1.15 

0.39 

Spd 

8.96 

10.23 

6.97 

5.18 

Num 

352 

224 

111 

17 
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Fig. 10. Same as Fig. 8 except for WVWV produced by Fig. 1 1. Same as Fig. 8 except for WVWV produced by 
CIMSS from the 7.3-jUm channel. CIMSS from GMS-5 over the western North Pacific region. 


did not evaluate WVWV; however, a later investiga- 
tion tested the sensitivity of GOES WVWV data on a 
spectral barotropic numerical hurricane track forecast 
model (Velden 1996). It was found that small improve- 
ments resulted in numerical track forecasts with the 
inclusion of the WVWV observations, despite the 
limitations of the analysis and forecast system employed. 

A natural extension to these previous track fore- 
cast studies is to assimilate the data into state-of-the- 
art, primitive-equation analysis and prediction 
systems in order to take full advantage of the satel- 
lite information and maximize the potential impact. 
Prediction systems developed at the Naval Research 
Laboratory (NRL) for operational use at the navy’s 
Fleet Numerical Meteorology and Oceanography 
Center (FNMOC) and from the Geophysical Fluid 
Dynamics Lab (GFDL) have shown significant skill 
in hurricane track forecasting. The NRL-FNMOC 
model (NOGAPS) is a global assimilation system that 


contains special tropical cyclone bogusing procedures 
(Goerss and Jeffries 1994). The GFDL model is a 
limited-area, multinested system that employs NCEP 
global fields as background information (Kurihara 
et al. 1995). Both systems were top performers in 
terms of objective hurricane track predictions during 
the prolific 1995 Atlantic season. 

These two prediction systems were employed to 
evaluate WVWV datasets prepared by CIMSS during 
the 1995 tropical cyclone season (Northern Hemi- 
sphere). Data were assimilated into NOGAPS for 
Atlantic Tropical Storm Chantal and western North 
Pacific Typhoons Ward and Angela. The Chantal 
datasets included winds derived at 1 2-h intervals from 
GOES-8 infrared (IR) and visible (VIS) imagery, and 
WVWV from the aforementioned three GOES-8 WV 
bands. The WVWV outnumber the IR-VIS winds by 
a factor of 6 and hence dominate the sample. The 
Ward and Angela GMS-5 WVWV datasets were pro- 
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cessed at 6-h intervals. The GFDL model was also run 
on the Chantal case and on Hurricane Gordon ( 1 994) 
using NCEP operational global analyses as the large- 
scale background fields, which were modified by the 
assimilated winds. The Gordon datasets contained 
only GOES-8 6.7 -fum WVWV, processed at 1 2-h in- 
tervals and were derived during the period when 
Gordon was upstream (south) of the dense U.S. upper- 
air observing system (Avila and Rappaport 1996). 

The results of the two-system prediction trials are 
given in Table 4. Track forecast errors were computed 
by using the observed track as verification and calcu- 
lating the distance between the observed and forecast 
locations. Mean forecast error improvement (ex- 
pressed as an overall percentage improvement over 
control runs that did not include the satellite data) and 
the number of improved forecasts relative to control 
runs are evaluated for each storm. Results are pre- 
sented for track predictions at 24, 48, and 72 h. 
Generally speaking, the data impacts are positive, in- 
cluding a significant reduction in forecast errors for 
the NOGAPS Chantal case (3 1 % at 48 h). An excep- 
tion is Supertyphoon Ward, which exhibited slight 
negative data impacts early on. In this case, the defi- 
nition of the broad upper-level outflow picked up by 
the WVWV was detrimental to the model steering. 
Overall, considering all cases from both model trials, 
assimilation of the wind information results in a re- 
duction of mean forecast errors by 6%- 1 3% over the 
24-72-h forecast period (all significant at the 95% 
level). The number of improved forecasts versus con- 
trol is also quite encouraging. 

It is informative to dissect the NOGAPS perfor- 
mance in the Chantal case since the most significant 
results were achieved there (all forecasts were im- 
proved at 48 and 72 h) and its datasets involved mul- 
tispectral winds. Table 5 presents the statistics for 
the control (CON) and satellite data enhanced (CON 
4- SAT) trials. The results are also compared to the 
climate-persistence (CLIPER) model. The reduction 
of mean forecast errors for the CON + SAT are con- 
siderable when compared to both the CON and 
CLIPER. The forecast error reduction is significant 
at > 92% level for all forecast periods. 

A satellite component analysis of the Chantal case 
is given in Table 6. In this set of trials, NOGAPS fore- 
casts were produced from analyses that excluded each 
individual component of the multispectral dataset in 
order to glean information on the impact of each data 
type. The results of the component analysis in this case 
reveal 1 ) the most significant contributor to error re- 


Table4. Numerical weather prediction data (WVWV) impact 
results for selected tropical cyclones. Mean track forecast errors 
(MFE) are evaluated and the results are expressed as a percentage 
improvement over the control (without WVWV) forecasts. GFDL 
is the Geophysical Fluids Dynamics Laboratory hurricane model 
and NOGAPS is the Navy Operational Global Atmospheric 
Prediction System. 


Forecast interval (h) 
24 48 72 

Hurricane Gordon (1994) 

GFDL 

15 

19 

5 

No. of cases (no. improved) 

6 (5) 

6(5) 

6(4) 

Supertyphoon Ward (1995) 

NOGAPS 

-6 

-4 

0 

No. of cases (no. improved) 

9(3) 

7(2) 

5(2) 

Supertyphoon Angela (1995) 

NOGAPS 

13 

6 

14 

No. of cases (no. improved) 

7 (5) 

7(6) 

7(4) 

Tropical Storm Chantal (1995) 

GFDL 

-9 

12 

12 

No. of cases (no. improved) 

10(4) 

8(6) 

6(5) 

Tropical Storm Chantal (1995) 

NOGAPS 

20 

31 

20 

No. of cases (no. improved) 

10(8) 

8(8) 

6 (6 ) 

Overall MFE 
improvement (%) 

6 

13 

10 

No. of cases (no. improved) 42 (25) 

36 (27) 

30 (21) 

Overall MFE (nm) 
for control 

94 

144 

197 

Overall MFE (nm) 
for control 4 sat 

88 

125 

177 


duction is the WVWV employed collectively, 2) the 
GOES sounder WVWV (7.0- and 7.3-jum observa- 
tions) is an important individual component at 24 h 
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Table 5. NOGAPS forecast impact results for Tropical Storm 
Chantal (1995). CLIPER is the climate persistence model, a 
standard for assessing skill. Mean forecast errors (MFE) are in 
nautical miles. 


Forecast interval (h) 
24 48 72 

Control (CON) MFE 

100 

179 

301 

Control -l- satwinds MFE 
(CON + SAT) 

80 

124 

242 

CLIPER MFE 

1 12 

210 

358 

% improvement 
(CON + SAT over CON) 

20 

31 

20 

% improvement 

(CON + SAT over CLIPER) 

29 

41 

33 

# forecasts 

10 

8 

6 

# improved forecasts 
(CON + SAT over CON) 

8 

8 

6 

T-test significance (%) 
(CON + SAT over CON) 

92 

97 

93 


and is also a significant contributor at 48 h, 3) the 
imager (6.7 fum) WVWV contributes significantly to 
the medium- and longer-range forecasts, 4) the high- 
density IR cloud-tracked winds are an important com- 
ponent at 24 h, and 5) the low-level visible winds do 
not appreciably contribute to the forecast improve- 
ment. This last point is likely a result of the fact that 
the VIS winds covered a very limited domain sur- 
rounding the storm that mainly just included the outer 
vortex region and were difficult to obtain in regions 
covered by the storm cirrus canopy. In addition, the 
NOGAPS contains a vortex bogusing scheme (Goerss 
and Jeffries 1994) that may have dampened the effects 
of these data. 

The revelations from the component analysis sup- 
port intuitive reasoning. The IR winds are most abun- 
dant in the storm region (cirrus canopy) and therefore 
contribute to short-term forecast improvements. The 
sounder WVWV provide information in the midlevels, 
closest to the storm steering current, and therefore 
have a considerable impact on the model track fore- 
casts. However, the limited viewing domain of the 
GOES-8 sounder confines the WVWV coverage and 


Table 6. NOGAPS sensitivity study of the impact of individual 
GOES-8 wind components in terms of the reduction on MFE for 
the Tropical Storm Chantal case. SNDRH20 is the combined 
sounder WVWV, IMGRH20 is the 6.1-jum WVWV, ALLH20 
is the combined imager and sounder WVWV, VIS is visible cloud- 
tracked winds, and IR is infrared cloud-tracked winds. Each 
component is subtracted from the collective dataset and the 
subsequent NOGAPS forecast results are shown. Mean forecast 
errors are in nautical miles. 


Forecast interval (h) 
24 48 72 


Control (CON) 

100 

179 

301 

Control 4* satwinds 

80 

124 

242 

Control + satwinds - SNDRH20 

94 

148 

246 

Control + satwinds - IMGRH20 

82 

151 

280 

Control + satwinds - ALLH20 

97 

156 

281 

Control + satwinds - VIS 

82 

133 

239 

Control + satwinds — IR 

90 

126 

236 


restricts the impact to the near environment in this 
case, hence, the shorter-term forecasts. The broad 
spatial coverage of the imager WVWV allows defi- 
nition of adjacent synoptic features and results in a 
contribution to forecast improvements at the medium 
and longer ranges (48-72 h). However, the primary 
message from this component analysis is that, in this 
particular case, the NOGAPS assimilation process 
benefits from the three-dimensional nature of the mul- 
tichannel observations, and the optimum results are 
achieved when the multispectral data are considered 
collectively. 

Additional model data impact experiments are be- 
ing performed at NRL. Also, recent changes in NCEP 
analysis and assimilation strategies (including the di- 
rect assimilation of satellite radiances) are leading to 
improved large-scale analyses. The ability of the 
WVWV observations to significantly contribute to the 
reduction in track forecast errors resulting from the 
GFDL model initialized with this new analysis sys- 
tem is being investigated. In the meantime, WVWV 
datasets are being produced in real time by CIMSS 
and routinely assimilated by FNMOC into the opera- 
tional NOGAPS. These datasets are also being used 
in qualitative analyses by the NCEP Tropical Predic- 
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tion Center in Miami and the U.S. military’s JTWC 
at Guam. 

7. Further applications 

The capability exists for current geostationary sat- 
ellites to provide nearly global WVWV coverage (ex- 
cluding polar regions). As discussed in the previous 
section, the operational utility of these data includes 
assimilation into global prediction systems, which 
should benefit from the increased information. The 
WVWV data can also contribute to research studies 
and operational problems involving regional analyses 
and local meteorological phenomena, particularly 
those occurring over marine areas. Applications to 
studies involving interannual scales are also possible. 
This section will briefly present a few examples and 
discuss other potential applications. 


a. Midlatitude oceanic analyses 

Forecasters have relied on subjective interpretation 
of water vapor imagery as an aid to analyzing meteo- 
rological conditions over oceanic regions (Weldon 
and Holmes 1991; Velden 1 987 ). The WVWV infor- 
mation can supplement this in two ways: quantifying 
the motion of features into actual wind speeds and as- 
sociating pressure heights with these features. Proper 
diagnosis of developing meteorological situations will 
benefit from this additional quantitative information. 

The eastern North Pacific is a region lacking in 
conventional observations, yet important to forecasts 
affecting U.S. interests. An example of WVWV cov- 
erage attainable from GOES-9 6.7-jUm observations 
is shown in Fig. 12. In this case, from 29 January 
1 996, a large evolving oceanic cyclone is present. The 
representation of upper-level flow around this system 
is well captured in the WVWV data. Case studies have 
shown that incorporation of these observations into 



Fig. 12. Example of WVWV from GOES-9 covering a strong extratropical cyclone in the eastern North Pacific on 29 January 
1996. Blue (yellow) vectors (kt) indicate winds with assigned heights above (below) 400 mb. 
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regional objective analyses can lead to improved 
medium-range forecasts (Velden et al. 1994; Velden 
1996). The potential addition of multilevel WVWV 
information from GOES-9 sounder bands over this 
region is particularly intriguing (Fig. 13) and is the 
focus of model impact studies under way. 

Improving the understanding of developing extra- 
tropical oceanic cyclones has been the goal of several 


meteorological field experiments over the past decade. 
The latest of these. Fronts and Atlantic Storm Tracks 
Experiment (FASTEX), will occur in 1997 (M. Shapiro 
1996, personal communication). It is planned that the 
intensive observing periods will include high-density 
WVWV fromGOES-S and Meteosat-5. Integration of 
these data with in situ aircraft observations, drop- 
sondes, ship rawinsondes, and satellite radiances will 
hopefully lead to an advance in 



Fig. 13. Example of GOES-9 WVWV derived from sounder bands at 7.0 pm (top) and 
7.3 ^m (bottom) on 19 April 1996. Yellow (beige) vectors (kt) indicate winds with assigned 


the scientific understanding 
necessary to enable detailed di- 
agnosis and prediction of the 
life cycles of North Atlantic 
storms. 

b . Tropical cyclones 

WVWV from GMS-5 are 
being routinely processed at 
CIMSS, delivered to the JTWC 
at Guam, and integrated into 
their subjective upper-level 
analyses. Figure 14 illustrates 
an example of the typical cov- 
erage achieved during this dem- 
onstrational support in 1995. 
Analysts at the JTWC indicate 
that the data are useful in depict- 
ing upper-level features and 
their evolution, which can play 
a role in tropical cyclone forma- 
tion and motion. Particularly in 
the vicinity of typhoons, the 
data depict upper-level structure 
and features not captured in the 
corresponding global analysis. 
This is illustrated in Fig. 15, 
which presents the WVWV for 
Typhoons Zack and Angela and 
the corresponding model guess 
winds (from the 6-h forecast of 
the NCEP operational global 
model, interpolated to WVWV 
location and pressure). This is 
an area lacking in conventional 
upper- air observations. Clearly, 
there are regions in the outflow 
pattern that are represented very 
differently. The impact of the 
WVWV winds on numerical 
forecasts of Angela’s track was 


heights above (below) 400 mb. The winds are plotted over their respective channels. presented in the previous section. 
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Fig. 14. Example of WVWV routinely produced for JTWC analyses from GMS-5 during Typhoons Zack (125°E) and Angela 
( 1 37°E) in the western North Pacific on 28 October 1995. Yellow (beige) vectors (kt) indicate winds with assigned heights above 
(below) 400 mb. 


The 1995 Atlantic hurricane season offered many 
opportunities to assess the GOES-8 WVWV product. 
Tropical Storm Pablo was a good example of quali- 
tative applications (see earlier discussion). Hurricane 
Erin presented another fine example. On 3 1 July, Erin 
was rapidly developing to hurricane status over the 
southeastern Bahama Islands. The official track fore- 
cast indicated a steady west-northwest movement with 
the likelihood of landfall over southern Florida within 
48 h. The actual track was characterized by a brief jog 
to the north-northwest on 1 August, which resulted in 
a landfall over central Florida near Vero Beach. The 
reason for this track perturbation was evident in 
GOES-8 water vapor observations. Image sequences 
indicated a subtle circulation to the northwest of Erin, 
traveling to the west. Implications are that Erin and 
this circulation became briefly involved in a binary- 
vortex interaction. This midlevel environmental circu- 
lation, while subtle, is captured in multispectra] WVWV 
on 1 August (Fig. 16). Operationally available data 


was not sufficient to describe this feature, resulting 
in analysis deficiencies and numerical guidance that 
did not indicate the jog in the track to the northwest. 

In addition to operational considerations, the WVWV 
fields are contributing to research analyses and studies 
aimed at a better understanding of the processes that 
can affect tropical cyclone motion and intensity. For 
example, the NOAA Hurricane Research Division 
(HRD) began synoptic flow experiments in 1982 that 
involve special aircraft reconnaissance missions to mea- 
sure storm environmental parameters (Franklin et al. 
1 996). These missions collect information from omega 
dropwindsondes below 350 mb (typical flight level). 
To supplement these datasets above 350 mb, CIMSS 
provides HRD with high-density multispectral GOES 
wind observations. From these datasets, objective 
analyses are being used in diagnostic case studies to 
examine the influence of synoptic-scale fields on 
changes in vortex intensity and track, with a goal to- 
ward improving short-term track forecasting. 
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Fig. 15. Plotted WVWV (yellow) and collocated NCEP Aviation model 6-h forecast wind vectors (red) during Typhoons Zack 
and Angela on 28 October 1995. Only vectors (kt) above 300 mb are shown. Note the differences in the typhoon outflow regimes. 


Tropical cyclone environmental steering is com- 
monly associated with a midtropospheric layer-mean 
flow (Pike 1985; Dong and Neumann 1986; Velden 
and Leslie 1991). Operational satellite-derived cloud- 
drift winds lack retrievable information in this part of 
the troposphere. In cloud-free regimes, WVWV can 
provide midlevel information, particularly if the 7.0- 
and 7.3-fjm observations are utilized. It was shown in 
the Tropical Storm Chantal example (previous sec- 
tion) that the ability to retrieve winds in the 400- 
600-mb layer can be an important contributor to 
defining the steering flow. 

It has long been recognized that the vertical wind 
shear of the horizontal wind can affect the structure 
and intensity of tropical cyclones (Gray 1968). More 
recent studies (Jones 1995; DeMaria 1995) have ex- 
amined shear in the context that the potential vortic- 
ity (PV) associated with the vortex is advected at 
different rates at different vertical levels. Tilting of PV 
in the vertical induces circulations that tend to restore 
the vortex to a vertical orientation and hence main- 


tain intensity. The ability of the vortex to survive ver- 
tical wind shear will depend on the Rossby penetra- 
tion depth (Hoskins et al. 1985) and the magnitude and 
profile of the vertical shear. DeMaria and Kaplan 
(1994) reported on an attempt to quantify vertical 
shear into a statistical hurricane intensity prediction 
scheme that has resulted in improvements over clima- 
tology persistence. However, the forecasts still explain 
only about 50% of the variability of the observed in- 
tensity changes. The vertical shear predictor in their 
scheme is the magnitude of the difference between the 
200- and 850-mb wind analyses. These two levels 
were chosen because most operational satellite- 
derived winds are located around these levels. The 
scheme may benefit from the multispectral WVWV, 
which will allow the evaluation of multiple levels. 

A real challenge to forecasters is the interaction of 
tropical cyclones with upper-level troughs, which can 
frequently affect both the storm motion and intensity. 
Water vapor image sequences have been used to iden- 
tify and diagnose these situations (Velden 1987). 
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Fig. 16. Plot of multispectral WVWV (kt) derived over marine areas from GOES-8 during Hurricane Erin on 1 August 1995. 
Yellow (beige) vectors indicate winds with assigned heights above (below) 300 mb. Note the circulation to the northwest of Erin that 
subsequently interacted with the tropical cyclone to alter its track. 


Traditional forecast philosophy suggests that these ployed toward model validation (Schmetz and 
situations purport to be hostile for tropical cyclones, Turpeinen 1988; van de Berg et al. 1991; Soden and 
given the presence of strong vertical wind shear com- Bretherton 1993; Schmetz and van de Berg 1994). 
monly found with the troughs. However, recent re- Significant differences have been noted between UTH 
search based within a PV framework suggests the fields derived from satellite observations and those in 
storm intensity will be modulated by the relative general circulation models. The natural extension of 
strengths and juxtaposition of the positive PV anomaly these studies is to utilize WVWV to infer upper-level 
associated with the trough, and the negative PV moisture transport and to diagnose seasonal or intra- 
anomaly associated with the outflow anticyclone annual changes linked to the UTH trends. 

(Molinari et al. 1 995). WVWV can play a key role in Schmetz et al. ( 1 995) describe the congruence of 
describing such interactions, and this is currently be- monthly mean fields of UTH and WVWV derived 
ing investigated. simultaneously from three geostationary satellites. 

The patterns (their Fig. 5) suggest a physical relation- 
c. Seasonal and climate scales ship in which the UTH fields are largely influenced 

Water vapor observations from geostationary sat- by upper-tropospheric transport. It is qualitatively 
ellites have been utilized to study short-term demonstrated that UTH and upper-tropospheric dy- 
climatologies. Specifically, the radiance information namics are directly observable on a nearly global scale 
from the 6-7.5-^im band can be transformed into an and that the UTH and WVWV fields are desirable 
upper-tropospheric relative humidity (UTH), which contributors to a climate-monitoring system, 
can then be related to atmospheric processes or em- Diagnosis of seasonal variations can also benefit 
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from the UTH and WVWV observations. For ex- 
ample, the 1994 Atlantic hurricane season was 
extremely inactive. Although easterly wave distur- 
bances emerging off the western African coast were 
numerous, only one evolved into a minimal hurricane. 
Certainly, cooler than normal SSTs in the region 
played a role, but it is likely that upper-level condi- 
tions also contributed to the lack of vertical development 
and organization of convection associated with the dis- 
turbances. To investigate the mean upper-tropospheric 
conditions present during the 1994 Atlantic hurricane 
season, UTH and WVWV fields were produced from 
Meteosat-3 observations nominally twice per day dur- 
ing the months of July-September 1994 and averaged 
over the 3-month period. 

The UTH fields were derived from Meteosat-3 
6.5-jum radiance measurements and a concurrent 
analysis of the state of the atmosphere (NCEP global 
analysis). The procedure utilized a forward radiative 
transfer model to compute radiances based on in- 
formation from the NCEP analy- 
sis. The observed radiances 
were compared with the com- 
puted radiances and a log-linear 
interpolation determined the fi- 
nal UTH value (Schmetz and 
Turpeinen 1988). The WVWV 
were also derived from the 
Meteosat-3 6.5-fim channel in a 
procedure similar to that de- 
scribed in section 3. From the 
3-month mean WVWV field, 
a layer-mean divergence was 
calculated. 

The respective mean UTH 
and divergence fields are shown 
in Fig. 17. The relatively dry 
upper troposphere between 1 2° 
and 30°N, stretching from the 
African coast to 65°W, is note- 
worthy. The layer-mean diver- 
gence analysis from the WVWV 
field indicates this region is 
characterized by upper-level 
convergence. Considered to- 
gether, these fields suggest 
mean subsidence was present 
over the region during this pe- 
riod. Such upper-tropospheric 
conditions may have acted to 
hinder convection associated 


with the emerging disturbances, thereby inhibiting 
further development. This analysis provides support- 
ing evidence that mean conditions over the eastern 
Atlantic were not favorable for development of west- 
ward-traveling disturbances emerging from Africa, 
which are a primary feature of the general circulation 
that can provide the initial stimulus for Atlantic hur- 
ricane development. A similar study involving analy- 
ses from the 1995 Atlantic hurricane season is under 
way in order to contrast conditions during this pro- 
lific hurricane season. 


8. Discussion 

Improved numerical weather analysis and predic- 
tion systems are placing increasing demands on ob- 
servation density and accuracy. Optimized data 
processing strategies must be developed in combina- 
tion with increased measurement precision in order 



Fig. 17. Composite of 3-month-mean (July-September 1994) upper-tropospheric relative 
humidity (color code in %) and divergence (contours, 10 -6 s -1 ) derived from Meteosat-3 
water vapor observations. Solid contours represent convergence. The divergence field was 
computed from a 3-month-mean WVWV field. 
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to meet these requirements. The remote sensing ca- 
pabilities from geostationary orbit are improving 
(Menzel and Purdom 1994). In this article, it was 
shown that initial attempts to optimize the data pro- 
cessing designed to take full advantage of these ca- 
pabilities can result in spatially coherent and useful 
measurements of atmospheric motion. In particular, 
wind vectors derived from water vapor image se- 
quences can complement our existing observing sys- 
tem and provide positive impact on both qualitative 
and quantitative analyses of the atmospheric state. 
Further research is necessary to refine the methodolo- 
gies (e.g., vector height assignments) in order to ad- 
vance the product quality and analysis impact. 

The NWP impact experiments presented in this 
paper and in previous studies, while limited, are en- 
couraging. Further studies are ongoing to assess 
the potential impact of the data on improving objec- 
tive analyses. The WVWV datasets can be reliably 
produced via automated methodologies at densities 
to meet varying meteorological requirements, in a 
time frame commensurate with operational needs. 
For example, large-scale (full disk) WVWV fields can 
be processed at reasonable densities for assimilation 
into global models. These fields can then be aug- 
mented with a higher-density product in meteoro- 
logically active regions (e.g., tropical cyclones) with 
an attendant quality control tuned to these active 
situations. WVWV datasets currently being produced 
at CIMSS from both GOES and GMS are being in- 
corporated into the navy’s operational global as- 
similation and prediction system (NOGAPS). At the 
time of this writing, such an approach using GOES 
WVWV was also being developed and tested by 
the NESDIS research division for NCEP operational 
considerations. 

A challenge to the NWP community is to continue 
to evaluate and explore ways to effectively assimi- 
late this data source. Traditionally, three-dimensional 
assimilation schemes have difficulty in retaining 
patchy, single-level information (e.g., cloud-drift 
winds, aircraft reports). Utilized collectively with a 
similar processing strategy from cloud-tracked winds 
(Velden et al. 1992; Nieman et al. 1997), the high- 
density multispectral WVWV should help retain sat- 
ellite information in the assimilation process by 
providing spatially coherent, three-dimensional cov- 
erage in the upper troposphere. Clearly, further study 
is necessary to explore the optimal assimilation of 
these data into objective analyses (i.e., observational 
weighting, balancing of the wind information with the 


mass fields, vertical blending, etc.). The ability 
of geostationary satellite sensing to provide high 
temporal interval datasets will also become an advan- 
tage to evolving assimilation strategies (e.g., four- 
dimensional variational schemes). 

As illustrated in this article, the WVWV fields can- 
not only contribute to operational objectives, but can 
benefit atmospheric research initiatives covering a 
wide range of scales. Understanding of meteorologi- 
cal processes will benefit as the data are integrated 
with the existing upper-air observational network and 
planned enhancements (e.g., autonomous aircraft, 
high-level jet aircraft with dropsonde capability, etc.). 

In summary, spatially coherent, high-resolution 
coverage of upper-tropospheric winds are possible 
from geostationary satellite remote sensing of water 
vapor. WVWV are comparable in quality to opera- 
tional upper-level cloud-motion vectors. However, 
vectors retrieved in regions lacking in definitive mois- 
ture structure are slightly less reliable, and the ability 
to obtain trackable information in extremely dry air 
masses and regions of strong subsidence is limited. 
The individual vectors represent single-level reports; 
however, upper-level vertical wind profiles are pos- 
sible if multispectral observations are employed (cur- 
rently possible only from GOES). The WVWV 
applications reported here have demonstrated the po- 
tential for analysis impact and make a strong case for 
enhanced remote sensing capabilities such as the La- 
ser Atmospheric Wind Sounder (Baker et al. 1995), 
and the Geostationary High-Resolution Interferometer 
Sounder (Smith et al. 1990) that are currently under 
consideration. Future improvements to the GOES 
system will include additional WV channels in the 7— 
7.5-jUm region on the imager, with improved horizon- 
tal resolution and radiometric precision. 

The next Indian satellite INSAT (to be deployed 
in 1997 or 1998) will join the geostationary meteoro- 
logical satellite fleet containing a water vapor imager. 
With the exception of polar regions, global high- 
density WVWV will then be possible. The interna- 
tional standardization and optimization of operational 
processing methods to enable coherent intersatellite 
data fusion is an objective of the World Meteorologi- 
cal Organization. 
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